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Abstract 

A recently proposed scaling law for the decay time of alpha particles is 
generalized for cluster decay. It is shown that for the decay of even-even 
parents, logTi/ 2 depends linearly on the scaling variable S=(Z c Zrf)°- 6 /- v /C7^ 
and on the square root of the reduced mass of cluster and daughter. 

PACS numbers: 23.90. +w,25.85.Ca,23.60.+e 
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In the last years alpha decay systematics have been theoretically reinvestigated |T|]- || 
partly as a result of the special interest in the studies of superheavy nuclei and cluster decay 
In the same time the cluster decay data has been accumulated J7| with a sufficiently high 
precision so that their systematics can be now studied almost model independently with 
guidance to the alpha decay. All the predictions or analysis of the cluster decay data given 
up to now 0]- jl6| are based on models. The prescriptions of these models and/or the 



physical meaning of the parameters used by them have been subjected to criticism. In this 
letter we present a model independent description of the accumulated data pointing out the 
most important variables which scale the cluster decay probabilities. 

The earliest phenomenological connection between the half lives and the Q-values of the 
alpha decays of radioactive series was proposed by Geiger and Nuttall |17| 

logT 1/2 = aQ- 1 ' 2 + b , (1) 

which proved to be very useful for the prediction of alpha decay half lives. The rule of 
Geiger and Nuttall has the disadvantage that the a and b parameters are dependent on the 
isotope chain. A universal scaling law for alpha decay half lives of the even-even parents has 
recently been proposed U 

logT 1/2 = 9.54^1 - 51.37 , (2) 

where is the charge of the daughter. It was shown that the known || alpha decay logTi/ 2 
of the even-even nuclei with Z>76 stay on this universal line with a rms deviation of 0.33. 
It is of interest to know if similar scaling law(s) exist (s) for cluster decay: 

A ?Z P ^ A ° < Z d + A c Z c (3) 

where the subscripts p, d, c refer to the parent, daughter and cluster, respectively. We look 
to even-even parents and clusters with the hope that the structure effects ]18j are limited 
and only the collective dynamics dominates the process. Relation (2) has been compared 
to what one obtains from that obtained from classical Gamow formula for the alpha decay 
constant 



A = uP P , (4) 

where us is the frequency with which the alpha particle exists at the barrier, P is the 
preformation probability (assumed to be constant in Ref. ]]]) and P is the barrier penetration 
factor assuming a square well plus a Coulomb potential for the radial dynamics. One obtains 

ZdZ r 



logTi/2 = C + 2 ■ log2 ■ d c e 2 yj2fx/h 2 [arccos (x) — x\J\ — x 2 } , (5) 

V Qa 



where 



Qc(R c + Rd) ,,,, 



and where Q c is the cluster decay Q-value and the R are the equivalent hard-sphere charge 
radii and C Q = log(ln2/ujP ). 

The scaling law (2) is not obvious from Eq. (5), but for x < 0.8, logP behaves ap- 
proximately linearly as a function of Z^ 6 / \/Q a . This is in fact the region of interest for 
the alpha decay of heavy nuclei. The above analysis indicates that the scaling variable is 
(Z a Z^) 6 / \JQ a . As a consequence, we have looked for the behavior of the known experimen- 
tal data on cluster decay (see e.g Table 1 from Ref. ||) as a function of the cluster scaling 
variable 

The data are presented in Fig. 1. There are only 3 known "chains" of cluster decay with 
more than one element: 14 C (logT^=11.02 from 222 Ra @, ^0^^=15.9 from 224 Ra @ and 
logT^=21.33 from 226 Ra 0), 24 Ne (logT^=20Al from 232 U 0, logT^=2AM from 
230 Th @ and logT^=25.2A from 234 U @) and 28 Mg {logT^=21.Q8 from 236 Pu and 



/o5'T 1 6 J|'=25.75 from 234 U [13|]). The experimental data concerning 14 C and 24 Ne cluster 
decay clearly shows a linear dependence as function of the scaling variable S. This analysis 
indicates a scaling law for the cluster decay (alpha included) similar to Eq. (2) 

logT 1/2 = C 1 (S-7) + C 2 . (8) 



The constant " seven" is subtracted from S in this equation simply so that the parameter C2 
has a numerical value which is close the actual experimental values shown in Fig. 1. The 
coefficients Ci and C2 can be extracted from the fit of the experimental data. The values 
of the Ci parameters are 6.3 for 4 He, 17.4 for 14 C, 20.7 for 24 Ne and 27.1 for 28 Mg. The 
corresponding values of C 2 are: -7.3, 8.0, 19.1 and 21.7. 

It is interesting to examine if there are some correlations between the Ci and C2 coeffi- 
cients corresponding to different cluster decays. Eq. (5) suggests an additional dependence 
on the masses. In the alpha decay case the dependence on y/Ji (the reduced mass) is very 
small: 0.5 % for the mass of the daughter (A^) in the range 150 to 250. For clusters heavier 
than alpha the dependence on yj~\x is very important. The analysis of the heavy cluster decay 
case is complicated by the fact that the preformation probability (prescission probability [|J ) 
plays a more important role in the majority of the theoretical models ||18|| - ||. There are 
also models |]7|JT6| which assure the prescission probability to be unity. 

Guided by Eq. (5) and by the fact that some models fll5],Pl indicate a JJl dependence of 
the prescission part also, we have plotted in Fig. 2 the Ci and C2 coefficients as a function 
of 



For the plotting purpose only we have used = 208 neglecting the small variation due 
to the different daughter masses. One can clearly see a linear dependence on y/JI of these 
coefficients (the fitted lines are given by C\ = 6.3y/JI — 6.2, C2 = 9.8y/Jl — 26.9). 

An alternative fit of the cluster decay data, in the spirit of those models [|7|,|T6l which 
consider this preformation probability equal to unity, can be performed with the help of Eq. 
(5). We have taken R c = 0.0354t4 c + 2.008 (fm) for the cluster radius which empirically 

1 /3 

reproduces the charge radii of the light clusters, and Rd = r A d ' for the daughter radius. 
The experimental data were fitted by using the two parameters C and r D . A 0.64 rms 
deviation from the experimental values has been obtained and the following values for the 
fitted parameters: C Q = —23.1, r D = 0.976 fm. The r D value is 0.25 fm smaller than the 



typical values for heavy nuclei, and this make the touching radius 1 - 1.5 fm smaller. This 
is an unreasonable reduction, but the "extra-penetrability" could simulate the preformation 
probability (the -23.1 value is consistent with a preformation probability P of about unity 
contributing to C ). 

Assuming that the preformation probability is different from unity, our analysis indicates 
that not only the postscission but the prescission dynamics also is dominated by the square 
root of the reduced mass. It is interesting to compare this conclusion with the prescriptions 
presented by different models. The present analysis is in accord with the prescription of 
Ref. and Q. The work of Blendowske and Walliser || indicates a linear dependence 
on A c (the cluster mass) of the prescission probability (the spectroscopic factor in Ref. ||), 
differing from the present conclusion. Barranco, Broglia and Bertsch |6j have obtained in 
their superffuid tunneling model a dependence of the prescission probability on the number 
of steps to the scission (an extra dependence of the gap parameter entering their formula on 
the mass of the parent nucleus does not affect this analysis). This number is very close to 
the reduced mass /i and not to ^//J, again differing from the our findings. 

The above analysis indicates a model independent law for the whole body of cluster decay 
data of the following form: 



done. Besides the 8 "in chain" cluster data considered in Fig. 1, 3 "single" cluster data have 
been taken into the fit: 20 O from 228 Th {logT^=20.9 0), 32 Si from 238 Pu (logT^=25.3 
H) and 34 Si from 242 Cm (logT^=23.2 @). The fit result gives ai = 9.1, h = -10.2, 
a 2 = 7.39, b 2 = —23.2, x = 0.416 and y = 0.613 with a 0.34 rms deviation of logTi/ 2 . 
Considering the important parameter x, a range of values from 0.4 to 0.6 can be obtained 
depending upon the various subsets of data used in the fit. A 0.58 rms is extracted for 
the heavy clusters only, which represents a fairly good description of the data if one has 
in mind that the largest deviation comes from 34 Si (/ogTi/2=24.45 as compared with the 




(10) 



A fit of the 119 alpha decays [0 and 11 cluster decays || from even-even parents has been 
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23.2 experimental value). This may indicate that the extrapolation of Eq. (10) to heavier 
clusters must be taken with caution. One would also like to see independent confirmation of 
the 34 Si experimental investigations [ 22"| . The apparent breakdown of this scaling law when 
going from cluster decay to fission could be understood by the fact that in the latter case 
the dynamics is not dominated by the Coulomb potential but by the collective potential up 
to the scission point. We compared our formula with model dependent results for heavier 
cluster decay like 48 Ca from 256 No. Our result is logT\j 2 = 27.9 while the result from Ref. 
19| is significantly smaller (logT\/ 2 = 18.9). Experimental information in this mass range 



are crucial. 

The scaling law, Eq.(10), can be straightforwardly used to produce tables with cluster 
decay half live predictions similar with those in Ref. [[HJ. Input parameters are the mass 
and charge numbers of cluster and daughter and the Q-value of the reaction. A detailed 
search through all the possible decays of the parents with 82 < Z p < 106 and clusters with 
2 < Z c < 20 shows the possibility to obtain experimental data for the decay of new clusters 
in this region: e.g 12 C from 220 Ra (logT 1/2 = 10.4) and from 222 Th (logT 1/2 = 10.08), 18 
from 226 Th (logTi/ 2 = 17.75), etc. To select these cases we have used similar constraints as 
in Ref. |T9"|| , namely logTi/ 2 < 28 and logTi/ 2 — logTi/ 2 (a) < 18. Only those nuclides for 



which the experimental masses are known |23| have been used. 



One can try to test the y/JI behavior in decays for which the daughter is different from the 
208 Pb region. The neutron deficient A « 120 region is particularly interesting. For example, 
the decay of 118 Ba into a 12 C cluster and 106 Sn. Eqs. (10) together with an experimentally 
extracted Q-value, Q c =15.10 MeV [pij, gives logTi/ 2 = 18.0. This extrapolation of the 



scaling law (8) to light parents gives significantly lower half lives as compared with other 
model dependent treatments (e.g. /ogT 1 / 2 =21.3 in Ref. [|ljj). Further experimental tests are 
required to validate one of these approaches. 

Eq. (10) represents the first model independent description of all known cluster decay 
data. The parameters a%, b\, a 2 , b 2 , x and y contain information on the dynamics of the 
decay. The actual theoretical models describing the cluster decay data are rather crude. 



Often their parameters loose their physical meaning, as for the unphysically small r Q dis- 
cussed above or e.g. the use of a zero-point motion energy even in the asymptotic region 
fl9|| . In our approach we have emphasized the most important variables (S, y//!) scaling 
the experimental data. We expect this new approach to be an important step toward a 
theoretical description of the cluster decay. 

In conclusion, we have obtained a new scaling law for the alpha and cluster decay of the 
even-even heavy nuclei. The scope of these scaling laws is to describe the regularities of the 
data, to put in evidence the peculiar behavior with respect to these regularities, to reveal 
the most important parameters entering the theoretical models and to guide new prediction 
for the cluster decay. New experimental data are necessary to further support the present 
analysis. 

MH and BAB would like to acknowledge support from the Alexander von Humboldt 
Foundation and NSF grant 90-17077. 
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Figure captions 



Figure 1 Experimental data for logTi/ 2 (sec) of the cluster decay of even-even parents 
as functions of the scaling variable S (Eq. (7)). The line noted by 4 He is given by Eq. (2) 
which represents the best linear fit to the experimental data. Other lines are drawn to guide 
the eye. 

Figure 2 Ci and C 2 coefficients entering Eq. (^) as function of ^JJL defined in Eq. (9). 
Lines represent the best linear fit. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9403008vl 




Figure 1. 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9403008vl 




Figure 2. 



